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Random hyperfine fields are essential to mechanisms of low-field magnetoresistance in organic semicon- 
ductors. Recent experiments have shown that another type of random field — fringe fields due to a nearby 
ferromagnet — can also dramatically affect the magnetoresistance. A theoretical analysis of the effect of these 
fringe fields is challenging, as the fringe field magnitudes and their correlation lengths are orders of magnitude 
larger than that of the hyperfine couplings. We extend a recent theory of organic magnetoresistance to calculate 
the magnetoresistance with both hyperfine and fringe fields present. This theory describes several key features 
of the experimental fringe-field magnetoresistance, including the applied fields where the magnetoresistance 
reaches extrema, the applied field range of large magnetoresistance effects from the fringe fields, and the sign 
of the effect. 



Introduction - A major thrust of research in organic spin- 
tronics concerns transport through organic semiconductors 
sandwiched between either magnetic^ ^ or non-magnetic^^ 
electrodes. Both classes of devices display large magnetore- 
sistance (MR), and have spurred significant experimentaP'S 
and theoretical^ interest. A complete description of the 
physics involved in this organic magnetoresistance (OMAR) 
is still evolving^^. For the non-magnetic case, one explanation 
of the MR relies on the existence of random, uncorrected hy- 
perfine fields (HFs) that vary from one localizing center to 
another (referred to as HF-OMAR).^ A recent experiment^^ 
has tested the importance and influence of randomizing fields 
by situating a single electrically isolated ferromagnet some 
specified distance from the organic semiconductor. When 
the ferromagnet is magnetically unsaturated, spatially vary- 
ing fringe fields (FFs) emanate outside the ferromagnet, and 
these FFs dramatically alter the MR lineshape (referred to as 
FF-OMAR). Such FFs could also be important in experiments 
with magnetic electrodes, such as organic spin valves A 
theoretical explanation of the FF influence has proven diffi- 
cult since the spatial distribution of FF's differs greatly from 
that of HF's; the magnitude and correlation lengths are both 
more than one order of magnitude larger. 

In this Rapid Communicatio n, we de monstrate that a re- 
cently proposed theory of OMARp^ESini^ based on percolation 
theory, offers a solution to the anomalous OMAR observed^ 
in the presence of FFs. The main ingredient is that field- 
dependent spin transitions (due to FFs and HFs) open other- 
wise blocked hopping pathways which in turn alter the sam- 
ple's resistance. We find that a careful study of the statistics of 
the FFs identifies several regions of behavior, including a re- 
gion close to the ferromagnetic film where FF gradients gov- 
ern the MR line shape, and a region farther away where HF 
re-emerge as an important element and where the dominant 
effect of the ferromagnet is through the FF magnitude. For 
regions closer to the ferromagnetic film, where FF gradients 
are important, we predict the size of the FF effect will not be 
very sensitive to the thickness of the ferromagnetic layer; as 
both the lateral size of the domains (and thus the field corre- 
lation length) and the magnetic moments will increase with 



thickness, the gradient should remain nearly constant. 




FIG. 1. (Color online) (a,b) Two examples of the ferromagnet 
with different domain configurations and normalized magnetizations 
(M/Ms) obtained by X-ray microscopy imaging in Ref. 14 (c) car- 
toon picture of two occupied sites and their local fields. Insets: exper- 
imental set-up and M/Ms versus applied field ^q. (d) Fringe fields, 
Bff^ and Bff^, in sl 2 jjm x 2 jjm area calculated 15 nm above the 
ferromagnet. 

We now give a short description and summary of the ex- 
periments on FF-OMAR reported in Ref. |T4| The or- 
ganic semiconductor device, called a semi-spin valve (Fig- 
ure [T](c)), consists of a ferromagnetic Co/Pt multilayered film 
with perpendicular magnetic anisotropy followed by a non- 
magnetic metal, a bottom electrode, an organic layer, and 
a Ca top electrode. FF-OMAR is observed even when the 
ferromagnetic layer is excluded from the current path, and 
therefore excludes effects such as spin-injection or tunnel- 
ing anisotropic magnetoresistance.^^ The reversal of magne- 
tization, M, in the ferromagnet occurs through nucleation, 
growth, and annihilation of magnetic domains. The satu- 
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ration magnetization, Ms, is 5.4 x 10^ A/m. The properties 
of the ferromagnetic films have been further characterized 
in Ref. 14 The organic layer is a 30- nm-thick film of 
tris(8-hydroxyquinolinealuminum (Alqs). The bottom elec- 
trode also serves as a spacer layer of variable thickness, Zi, to 
separate the magnetoresistive material, Alqs, from the mag- 
netic layer that is the source of the magnetic FFs. The FF's 
strength and spatial-correlation length depend on the magnetic 
domain configuration as well as on the distance from the fer- 
romagnet to the organic film. Typical experimental data is 
shown in Figure [2] (a). Unlike ordinary HF-OMAR, the FF- 
OMAR is hysteretic and extends to a much larger field-scale. 
In Ref. 14 transmission x-ray microscopy (TXM) based on the 
x-ray magnetic circular dichroism (XMCD) effect was used 
to determine the Co/Pt layer's microscopic magnetic domain 
structure as a function of the applied perpendicular field (see 
Figure [T](a,b)). We obtain here, through methods described 
below, the FF distribution at a given distance above the fer- 
romagnetic film from magnetostatic modeling (see Figure [T] 
(d)). 

Theory - In disordered organic films transport occurs by 
hopping between localizing sites. The fermionic nature of the 
charge carriers forbids the formation of doubly occupied sites 
in a triplet state. Hence a hopping polaron experiences a re- 
duction in the number of accessible sites which affects the 
transport properties of the organic semiconductor. This leads 
to the formation of bottlenecks where a transport pathway is 
restricted due to the spin-blocking described. The situation 
can be alleviated when two polarons on nearby sites (a po- 
laron pair) experience local magnetic fields which allow the 
charge-blocking triplet polaron pair to undergo a transition to 
the singlet state which allows for a doubly occupied site to 
form. However if the spin transition is slow enough, it may 
be more expedient for the blocked charge to bypass the other 
occupied site. This produces a competition between two pro- 
cesses: formation of a doubly occupied site (e.g. a bipolaron 
or a doubly occupied deep trap) and disassociation of the po- 
laron pair. For local fields which are due to the nuclear spins, 
theory predicts the cross-over between the two processes can 
be seen in the MR when the hopping rate, vq, is varied from 
slow hopping to fast hopping (compared to the strength of the 
HF).li^ 

Solution to the spin-independent transport problem in such 
systems comes from percolation theory. The percolation 
theory solution of OMAR is encapsulated in the fol lowing; 
equation for the threshold hopping distance, rc| ii|i6 | i7 | 

j\nNeffr^dr = Bc (1) 

with Neff =N — Nt^ apT^s^T where Nt = \N\ is the den- 
sity of polarons that form a triplet spin state with some specific 
polaron; N\ is the concentration of injected polarons with lo- 
calization length I. N is the density of sites in the system 
and considered to be M (i.e. dilute carrier concentration 
approximation), a is a number less than unity. Eq. ([T]) is spin- 
dependent through the probability for a triplet polaron pair 
to transition to a singlet polaron pair, pr^s- is the average 
number of sites within a distance of one another that exist in 



the percolating cluster; in three dimensions, 5^ ^ 2.7 as found 
from numerical simulations.^^ The theory is valid for unipo- 
lar transport, though many of its qualitative features can also 
be applied to the bipolar regim^^^^ probably relevant for the 
FF-OMAR experiments^^. The theory has been further bol- 
stered by success in explaining features of experimentJSl in 
the bipolar regime. 

After defining MR= {\Rc{B ) - Rc{^ )\ / Rc{G)) with Rc = 
RQe^'^l^, we use Eq. {I} to writJ^IISlIll 

-Y-j pill rTZ rye 

MR ^ ^ / / / {ps{B)-ps{^))y^^mMydM^, 
6ny^ Jo Jo Jo 

(2) 

where y = r/£ and we have made the substitution pr^s = 
5 [1 — Z^^"] in terms of the easier to evaluate quantity ps^s = Ps 
and r\ = Nj /N?^ Estimating r\ is difficult since the carrier 
concentration is not known with precision and deeply charged 
traps likely exist that do not contribute to the carrier popula- 
tion but can contribute to the MR.-^^^ It has been established 
thaP 

where Ps is the singlet projection operator, m and m' are in- 
dices denoting the eigenstates of the total Hamiltonian, CO rep- 
resents the eigenvalues, and Xh = ^ exp(2r) is the hopping 
time between two sites separated by a distance r. For con- 
creteness, we choose I = 0.2 nm and = 1 nm throughout 
this Rapid Communication.ESlll We assume vq = 10^ in fre- 
quency units of 1 mT. 

The system Hamiltonian of the polaron pair, = J^o + 
J^ff is composed of terms that are responsible for 'normal' 
OMAR, 

JTo = {(Hhfin ) + G)o) • 5i + (G^hfiri) + «)o) • 5*2, (4) 

and a new term arising from the FFs: 

J^ff = (Offin ) • 5i + G)//(r2) • ^2. (5) 

Both HFs and FFs have spatial dependence; however as shown 
below, the dependences are very different. HFs are uncorre- 
lated from site to site (on the order of 1 nm) and their mag- 
nitudes follow a Gaussian distribution of width a which is on 
the order of a few milliteslas. FFs are correlated over a larger 
distance (~ 100 nm) and their magnitudes are ~ 100 mT for 
a spacer length of Zi = 15 nm.^"^ 

Calculation - The existence of fringe fields adds consider- 
able complexity to a calculation of Eq. ^ since spatial vari- 
ables must be accounted for in the eigensystem of ps. Ad- 
ditionally the placement of a bottleneck could be in either a 
region of small or large FF in comparison to the HE; also the 
FF may either be constant or sharply vary between the two 
occupied sites. These issues are not encountered when the lo- 
cal fields are entirely uncorrected (such as was the case with 
HE). To account for these complexities, during the integration 
we sample pairs of sites (i.e. bottlenecks) randomly through- 
out the organic semiconductor volume. The spatial integrals 



3 



o 



Experiment of Wang et al. 
100 nm 



"^■50 nm^i^ 1 ^ 


. 20 niTiYM/ 













500 1000 -10 




1.0 
0.8 
0.6 
^ 0.4 
^ 0.2 
0.0 
-0.2 



Bo (mT) 


Bq (mT) 


(c) 


z/=100nm_ 




/ Z/=65 . 




/ y/^ i^45 nm , 




— ""^^l Zi=15 nm 



-50 



50 




).001 0.01 0.1 



10 100 1000 



Bq (mT) 



M, (5.4x10^ A/m) 



FIG. 2. (a) MC as measured for different spacer distances by Ref . "T? 
(b) HF-OMAR calculated from our theory, (c) FF-OMAR calculated 
from our theory for several spacer distances with a = 3 mT and d = 
10 nm. A moving average is taken to smooth the curves. The applied 
field range is limited by the availability of X-ray images at larger Bq. 
Thick arrow indicates the direction of field ramping. 



FIG. 3. (Color online) (a) The quantity P{Bo) (Eq. [6]) at 5o = mT 
(dotted) and 36 mT (solid) for three different hyperfine couplings 
a = 0, 1, and 4 mT. (b) The calculated % MR for the same hyper- 
fine couplings at 36 mT. The field 36 mT is chosen since it nearly 
corresponds to the magnetization M/Ms = 0. 



(as well as the six integrals resulting from the average over 
hyperfine configurations) are done numerically by the Monte 
Carlo integration method: 

r2n rn rye 

P{Bo)= / / {ps{y,e,(\>))y^smddyddd(^^ 
Jo Jo Jo 

1 ^ 

where the random numbers are chosen uniformly from the 
ranges < y^ < yc, < < tt, and < < 2n. Angu- 
lar brackets denote the averaging over hyperfine field distri- 
bution; the coordinates of the two hyperfine field vectors are 
also random variables included in the summation. Conver- 
gence is slow; for the results presented herein, ^ = 1 — 10 mil- 
lion. Whenever the applied field is changed, the calculation 
requires the FFs to be recalculated given the new domain con- 
figuration. There are 3 1 X-ray images of the domains for 3 1 
different applied fields within the magnetic switching regime. 
We calculate the MR at these applied fields in an 1 /nm x 1 /nm 
X d volume. Larger volumes account for the domain struc- 
ture more accurately but are computationally expensive. In 
the field-range where the ferromagnet is saturated (no FFs), 
we calculate MR from the HFs only. 

Results- Figure |2jc) shows our main results for the MR, 
to be compared with the experimental results from Ref. [T4l 
Figure |2ja). Note that Ref. 14 's results are for the magneto- 
conductivity (MC), which implies that their MR is of opposite 
sign to our calculations (when a MR is calculated for unipolar 
transport, the MR seen in bipolar transport usually has oppo- 
site sign due to the fact that forming a doubly occupied site. 



i.e. an exciton, reduces currenPSEHl)^ xhe experimental and 
theoretical curves exhibit similar trends as the distance be- 
tween the organic and ferromagnetic layer increases. For ex- 
ample, a region of negative MR is seen at small positive fields 
for small distances, and the area of this region decreases as 
the distance increases. For even larger distances eventually 
HF-OMAR should emerge as FFs become weaker. This is ar- 
tificially demonstrated in Figure |3](b) by controlling the mag- 
nitude of the FFs by adjusting Ms. As the FFs (as well as 
the size of their gradients) are reduced, HF-OMAR is retained 
(~ 18.5 % in this case). By adjusting the HF coupling, it be- 
comes apparent that there is a competition between FF and 
HFs in determining the overall behavior. The X-ray images 
exist for only ramping up the field so the hysteretic nature of 
the FF-OMAR cannot be seen from the calculation. 

A natural question arises regarding the minima observed in 
the MR line shape. Why and where do the minima occur? We 
find that this feature corresponds well to behavior observed in 
the FF statistics. From the analysis of the FF statistics per- 
formed in Ref. 14 it is known that the FF correlation length 
was two orders of magnitude larger than that of the HF cor- 
relation length (~ 100 nm versus 1 nm). However, possibly 
a better indication of the efficacy of FFs in altering the MR 
is the FF gradient between the bottleneck sites. We find that 
the gradients can be comparable to the HF gradients so it is 
reasonable to expect the FF to modulate the MR. 

Our statistical analysis is depicted in Figure]?] where the av- 
erages and standard deviations of FFs are determined. Several 
conclusions can be derived from the FF statistics. One key 
point is that even at = mT, there still exists large FFs 
such that Bff ^ B^f. This is shown in Figure [41(e) for the 
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FIG. 4. (Color online) (a) and (b) are two examples of the partial 
derivative of Bff^ with respect to z distributions, (b) The width nar- 
rows as distance from the ferromagnet increases, (b) The distribu- 
tions are wider near M = which corresponds to 5o ~ 36 mT. (c) 
An example of the total field (minus HF) distribution at Bq = mT 
for small and large z/. (d) Standard deviations, a^z, of d^bzir) com- 
puted at random points in the specified volume with bottom edge at 
Zi. Plots (e) and (f) follow the same direction of labels for different 
Zi. Oij for {ij} G {x,y,z}) display similar behavior. The average 
gradients (not shown), e.g. d^bzir), do not stray significantly from 
zero, (e) The average total field (minus HF) magnitude as a function 
of applied field for several different Zi. (f) The standard deviation of 
the total field (minus HF) magnitude as a function of applied field for 
several different Zi. The vector B is Bq +^//- 

smaller Zi where Bff ^ 100 mT at = mT (zi = 15 nm). 
While the width of the distribution is ^ 30 mT, the FFs still 
quench practically all the HFs. In view of Figure [4] (f) also, 
the situation is unchanged; the average total field increases as 
does the width but still all HFs are overpowered. This sug- 
gests that the mechanism that produces the deviation from the 
Lorentzian line shape — referred to here as 'ears' — is not 
due to a mixture of HF-OMAR components with uniform FFs 
adding to the applied field (this is the "uncorrected B -field 
model" of Ref. The 'ear' cannot be due to any decrease 
in net field; B increases while Gb decreases modestly (Figure 
|4](e,f)) and an increase in total field tends to result in MR > 
which is not observed in that field range. However if the dis- 
tance Zi increases, the situation changes because the average 
total field reduces dramatically, as demonstrated in Figure |4] 
(c). For such cases, the average field 5 in a small finite field 
can actually be larger than the average field at zero field; this 
effect reduces the resistance and shows up as MR < (not 
shown). 

The statistical analysis of the previous paragraph points to 



the OMAR 'ear' for low Zi being due to the gradients in the 
FFs. Figure |4] (a) shows that typical sizes of the FF gradi- 
ents can be up to a few millitesla per nanometer which would 
make the effect larger or at least comparable to HF-OMAR. 
The trend in a^z with applied field correlates with the 'ear' 
which lend further credence to the importance of the FF gra- 
dients. This is explainable from the fact that at M = 0, there 
are equal numbers of up and down domains which give rise to 
maximally varying FFs; in short, the probability that the two 
sites of the bottleneck feel local fields that are different is a 
maximum at M = 0. As the magnetization approaches satura- 
tion, one expects the standard deviation to be zero as the FFs 
are quenched - however the quality of the X-ray images di- 
minishes in such cases and the domain configurations cannot 
be extracted for M/Ms nearer to unity. HF-OMAR does not 
immediately emerge when the FF gradients get smaller than 
the HF gradients - 65 nm curve in Figure [4] (d) is very nar- 
row so FF gradients are small. However the MR line shape is 
still significantly modified from that of HF-OMAR (Figure [2] 
(b)). While the gradients are small, the HFs still are influenced 
by the FF (even at zero applied field), as the FFs, though not 
varying spatially much, are still larger than the HFs. 

Figure|3](a) shows the quantity P{Bq) for three different hy- 
perfine couplings, a. First consider non-zero HF and the = 
mT lines (dotted): for the smallest M^, the FFs are negligi- 
ble so P and MR are independent of Ms. When Ms increases 
to where B ^ a, the FFs still have very small gradients (e.g. 
the Czz of Figure [4] will scale linearly with so at = 0.01, 
Gzz ~ 0.03 mT/nm <^Ghf mT/nm) so the FFs effectively 
act as uniform fields across a site pair. This shows up as an 
increase in P{0) and a decrease in MR(36 mT) since even 
though the applied field is fixed, increasing Ms is effectively 
increasing the uniform field felt by a site pair. Eventually Ms 
can be increased enough to where Gzz becomes comparable to 
Ghf', the variations in FFs between the site pair can now cause 
spin transitions more effectively. Since the FFs random vari- 
ations are larger at 36 mT than at mT, the MR is negative. 
For larger a, the P{0) and MR curve shift to higher Ms since 
Ghf is larger. When a = mT, a different picture emerges; at 
Ms =0, there can be no spin-transitions between singlets and 
triplets so ps = I and P{0) — ^ 523.6. In the opposite 
limit of large Ms, ps = 1/2 so P{0) converges to ^ ^yl- 

Conclusions- The present approach has shown that the 
salient features of FF-OMAR are explainable by our recent 
theory of OMAR based on calculating magneto-transport with 
percolation theory. This work serves as a springboard for fur- 
ther calculations to be performed within other OMAR theories 
- especially those that handle ambipolar carriers.!^ To find bet- 
ter agreement with experiment requires studying the depen- 
dence of MR on the hopping rate, vq. That task necessitates 
more intensive numerical work and thus has been avoided in 
this current discussion. This line of inquiry is also expected 
to elucidate how the FF-mechanism is similar to either the HF 
or the Ag mechanisms. 
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